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Abstract 

A central aim of physics is to describe the dynamics of physical sys- 
tems. Schrodinger's equation does this for isolated quantum systems. 
Describing the time evolution of a quantum system that interacts with its 
environment, in its most general form, has proved to be difficult because 
the dynamics is dependent on the state of the environment and the corre- 
lations with it. For discrete processes, such as quantum gates or chemical 
reactions, quantum process tomography provides the complete descrip- 
tion of the dynamics, provided that the initial states of the system and 
the environment are independent of each other. However, many physical 
systems are correlated with the environment at the beginning of the ex- 
periment. Here, we give a prescription of quantum process tomography 
that yields the complete description of the dynamics of the system even 
when the initial correlations are present. Surprisingly, our method also 
gives quantitative expressions for the initial correlation. 

There is a rich history to the studies of decoherence of quantum systems due 
to the interactions with the surrounding degrees of freedom. When the dynamics 
of the system (S) is Markovian it can be described by a master equation [TJ|2J|3]. 
Nowadays many researchers are interested in systems that are non-Markovian, 
as there is mounting evidence that some natural systems of importance may 
be non-Markovian £Q and such features may allow to manipulate and control 
quantum systems in desired ways. There is also a great deal of interest in 
systems that are initially correlated with their environments (£) because non- 
Markovianity and initial system-environment (S£) correlations are intimately 
related [51 HE]. 



Grasping the mathematical and physical aspects of non-Markovian systems, 
especially with initial S£ correlations, has proved to be a tough road. Nev- 
ertheless, there is a great deal of progress on deciding whether a system is 
non-Markovian in the recent years [El El EH E] • However, avoiding the initial 
S£ correlations is not always possible in reality [TH [T31 E] • Working with initial 
correlations in practice has proved to be much trickier than in theory. This is 
because the presence of correlations do not allow for a clear definition of the 
state S independent from the state of £ and vice versa. Physical systems are 
complicated and have many additional degrees of freedom that are not of exper- 
imental interest. Yet these extra degrees of freedom interact with the degrees 
of interest leading to correlations. Therefore initially uncorrected S£ state is 
often an approximation. 

In theory of open quantum systems, discrete quantum transformations are 
described by the dynamical map formalism [HI [16]: B{p s ) = pf . The dynamical 
map can be thought of as coming from the contraction of S£ unitary dynamics. 
Let us write the state of S£ as 

p ss = p s ®p s + x S£ , (1) 

where \ S£ 1S the correlations matrix [17] . The dynamical map is the map- 
ping from the initial states of S to the final states of S, resulting from unitary 
dynamics of the S£ state 

B(p s ) = pf = tr £ [U P S£ U^] (2) 
= tr £ [Up s ® pFrf] + tr £ [U X S£ U^] (3) 
= B CF (p s )+B^, (4) 

where B CF is a completely positive map and B aS is the affine correction term due 
to the initial S£ correlations. This means that B may not a be completely pos- 
itive map when \ S£ 7^ 0, nevertheless it fully describes the dynamics of S [19]. 
However, to determine such a map experimentally would require preparing dif- 
ferent states of S while keeping the S£ correlations fixed. Such preparations are 
not operationally feasible because altering the state of S will also alter the S£ 
correlations. Therefore, a nonpositive dynamical map is not an operationally 
meaningful quantity. 

The operational approach to quantum dynamics relies on the fact that quan- 
tum theory is a theory of preparations and measurements. The experimental 
method to determine a dynamical map corresponding to a quantum process is 
called quantum process tomography (QPT) [201 HI]- It is the central tool in 
determining a discrete quantum process; e.g. quantum gates [2"2l [231 |2"H 1231 l2"6l 
[27l [28] [29l 130] or chemical reactions [311 E21 EH] • To see the difference between 
QPT and dynamical maps let us review the four basic steps necessary to carried 
out QPT [31 [35]: 

(i) Input states that span the space of S are prepared. 

(ii) The input states are sent through the process. 
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Figure 1: Standard 
quantum process tomography. At the beginning of the experiment the system- 
environment state is uncorrelated. A preparation (.4) is made on the system 
and the corresponding output state Q is observed. This process is described by 
the completely positive map of Eq. Q, which is a function of initial state of 
environment and the unitary dynamics. It maps the initial states of the system 
to output states Q. 



(iii) The corresponding output states are determined by quantum state tomog- 
raphy. 

(iv) The knowledge of input states, the corresponding output states, and as- 
suming linearity completely determines the process. 

Let us denote input states as P and output states as Q. The first step of 
QPT is state preparation. A preparation procedure takes an unknown state of S 
to a known state of S. Mathematically, it is described by a completely positive 
map acting on the system [3B]. For instance, consider a set of preparations that 
project S into pure states: [A^ m) 8 2] (p S£ ) = P^ <g> /K m ). Since P( m ) is 
a pure state, the post-preparation S£ state is fully uncorrelated, where p e \( m > 
is the conditional state of £ . I is the identity operator acting on £, as we 
assume that the preparation procedure only acts on S and not £ . We will 
discuss the implications of relaxing this assumption in Discussions. Lastly, if 
the preparation is not trace preserving, it should be divided tr [^"^(p 5 ^)] for 
normalisation. 

The S£ evolution, after the preparation yields the output state: 



Q (m) = trg 
= tr f 



U[A {m) ®l]{p S£ )uA (5) 
UP {m) ®p £ \ {m) uA . (6) 



The key difference between the dynamical map in Eqs. ([2| and ^ is the act of 
state preparation. Because dynamical maps of do accommodate state prepara- 
tion, they are not operationally defined. In the presence of initial S£ correla- 
tions, state preparation affects the state of £ in a nontrivial manner. That is, 
the state of £ in Eq. ^ is conditioned by the choice of the preparations. 
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In deriving the standard QPT procedure it is implicitly assumed that the 
initial state of S£ is uncorrelated [37], i.e., the state of £ is thought to be a 
constant of the problem. When that is the case, the state of £ in Eq. ([6| is 
not conditioned by the preparation procedure. In this case the derived map for 
the process is completely positive and is the same as the completely positive 
dynamical map in Eq. Q. See Fig. 1 for a graphical illustration. In the 
presence on initial S£ correlation, the conditional state of £ will be different 
for each preparation, and the assumption of linearity in step (iv) of QPT is 
violated, i.e. the map is a function of the preparation procedure. Such maps 
are nonpositivc, nonlinear, or simply put nonsensical [341 135] . 

It then begs the question, can we determine the dynamics of a system that is 
initially correlated with £ ? This is an important question for two reasons: First, 
there may be physical system of interest that may have initial correlations. Is it 
possible to study their dynamics? Second, for foundational reasons we may care 
to know what are the limitations in describing the dynamics of physical systems. 
A partial solution to these questions was given in [331 EE] • In this article we 
show that not only complete dynamics of initially correlated system can be 
determined, we can also determine the contribution due to initial correlations. 



Results 

A map on a map 

QPT is performed out by noting how input states, that span the space of S, 
map to output states. The key insight in what follows is that it is not the input 
states of S that are relevant, rather it is the preparation procedures itself, i.e., 
the preparation map A. For a dg dimensional system there are dg linearly inde- 
pendent states that span its space. However, there are d$ linearly independent 
operations (preparations) that span the space of preparations. If we determine 
the corresponding output states for a set of linearly independent preparations 
then by linearity we have can predict the output state for any preparation. Let 
us denote this map as A'f-map. 

The form of A^-map arises naturally when considering the whole process 
in physical terms: At the beginning of the experiment S£ is in an unknown 
(correlated) state, p . The system is prepared into a known input state by the 
preparation procedure A, followed by a joint unitary dynamics. The output is 
given by tracing over the environmental degrees of freedom: 

Q = tr £ [U[A®1]( P S£ )U^] . (7) 

We want a map acting on the preparation map A and yielding the output state 
Q: M {A) — > Q. Then Al-map is everything on the right hand side of Eq. Q 
that is not A. The expression for M-map in terms of matrix indices is 

M-rr'r";ss's" = U re y a pf"a,s" pU* es i p. (8) 
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Figure 2: Quantum 
process tomography with M-map. At the beginning of the experiment the 
system-environment state is correlated. A preparation is made on the system 
and the corresponding output state Q is observed. This process is described 
by the completely positive map A4, which is a function of the initial system- 
environment state and the unitary dynamics. The A4-map takes preparations 
A to output states Q. 



Above a sum over repeated indices is implied. A4-map is a 'super super-operator' 
that acts on the super operator A. AW-map is a d$ x <i| tensor, which is 
contracted with a preparation A, a d 2 s x d 2 s tensor, yielding the output state Q, 
a d$ x d$ matrix. In term of matrix indices, the action is as follows: 

Qrs — A^ Tr i T " - : ss's" (Ar'r":s's" ) • (9) 

Again, a sum over repeated indices is implied. In Methods, a full derivation for 
A4-map in the last equation is given. See Fig. 2 for a graphical illustration of 
AW-map. 

Note that, in standard quantum process tomography state of £ is a constant 
of the process, here it is the initial S£ state that is the constant of the process, 
i.e., it is a fixed quantity. Physically, the constancy of p S£ means that the 
experiment should be initialised in the same manner for every run, and then a 
preparation on S can be made. 

A4-map contains both U and p ; however knowing M. is not sufficient to 
determine U and p S£ . As expected, it should not be possible to determine U 
and p S£ through measurements and preparations on the system alone with- 
out access to the environment. Conversely, AW-map contains all information 
necessary to fully determine the output state for any preparation of S. The ad- 
vantage of dealing with the A^-map is that we have separated the preparation 
procedure from uncontrollable dynamical elements and the initial conditions. 
A'J-map contains all of the dynamical information for the system and in the 
next section we will extract some of this information from the A4-map. First 
let us mention some properties of A'f-map derived in Methods: Its action on 
a mixture of preparations is linear, it preserves trace, it preserves Hermiticity, 
and it is completely positive. 

In Methods wc show that AW— map can be experimentally determined by 
making a set of linearly independent preparation of the system. This is similar 
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to what one has to do in standard QPT. In standard QPT a linearly indepen- 
dent set of states are fed into the process and the corresponding outcomes are 
observed. Knowing the inputs and the outputs the standard process map is de- 
termined. The difference here is that a linearly independent set of preparations 
are fed in to the process. This is of our major result of this paper: We have 
given a prescription to determine the dynamics of a system in an operational 
way, i.e., a mapping from preparations to output states. 

Quantifying initial correlations 

The .M-map contains the dynamics of the system before any preparation is made 
on the system. It is a function of the initial state p state as well as the S£ 
unitary transformation. .M-map is a tensor, taking its trace with respect to the 
indices that belong to the initial state of S we can obtain the dynamics of the 
system as if the initial correlations we absent. Using this with the knowledge of 
the initial state of S, in Methods we show that from M. we can derive another 
matrix, 

£rr'r";ss's" = '/^ Ure,r'aPr" s" Pad^se.s' 0' 0-®) 
a/3e 

Matrix C is fully determinable from M.— map and the two are the same when 
there are no initial correlations. We will call the difference between A4 and £, 
K, = M — C, the memory matrix: 

k-rr'r";ss' s" = ^ ] Ure,r'aXr"a;s" U se,s' fS ■ (H) 
apt 

Since M. contains p S£ and C contains p s ® p £ , the difference between the two 
is a function of only x ■ The action of the correlation-memory matrix on a 
preparation yields 

X s A (t) = KL{A) = tx £ [U[A®l](x SS )U% (12) 

which is the coherence coming into the system from the initial correlations. 
For non-Markovian dynamics the future state of S may depend on the initial 
S£ correlations. This is the non-Markovian 'memory' due to the initial S£ 
correlations and it is a key feature of non-Markovian dynamics [7] . 

The correlation-memory matrix is an important result for studying non- 
Markovian systems. It is an operational way of measuring the information that 
flows into S due to correlations at the time of the preparation. Once M.— map 
is determined, we have the full knowledge of the dynamics of S that is due to 
the initial correlations. The correlation-memory matrix provides quantitative 
information about the initial correlation and it is more than a witness for initial 
correlations [12] . 
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Operational meaning of not-completely positive maps 



For the special case, when the preparation is chosen to be the identity map, we 
get pure dynamics of the correlation-memory matrix 

K{l)=tr £ [U X S£ U% (13) 

which is the reduced dynamics of S£ correlations. This is exactly B aS in Eq. |ZJ 
From M-map we can determine matrices C and IC. In turn, from C we can get 



B CP (see Eq. (28 1) and from JC we can get B aS , and together they give us B of 
Eq. which can be a not-completely positive map. This gives not-completely 
positive maps an operational meaning. 



Discussion 

.M-map is the result of a quantum process tomography procedure for initially 
correlated system-environment states. It is acts on the preparation of the ini- 
tial state of the system, and only contains dynamical information. We study 
the properties of .A/f-map, showing it to be linear, preserving of trace and Her- 
miticity, and completely positive. Dynamical information about the evolution 
of the initial correlations can be retrieved from M— map, in the form of the 
correlation-memory matrix IC. .M-map allows us to determine the output state 
for any preparation of the system, while the correlation-memory matrix IC pro- 
vides a quantitative expression for the coherence due to the initial correlations. 

An important question is when is .M-map relevant? Clearly, when S and £ 
are initially uncorrelated then IC will be zero. Alternatively, just the presence 
of initial S£ correlations does not warrant for VW-map. Suppose \ S£ but 
trc\Ux S£ U^] = 0, then the completely positive map of Eq. Q would suffice to 
describe the dynamics correctly for any preparation of S |18| . 

One downside to A4~ map is that it requires a lot of resources to construct. 
In standard quantum process tomography d s input states are fed through the 
process and the corresponding output states are determined. To determine 
M. map, d s preparations are necessary, which is a significant growth over the 
standard procedure. Therefore an efficient way, such as compressed sensing [391 
140) . to determine this map is desirable. This should be possible, as determining 
A'J-map is equivalent to carrying out ds standard quantum process tomography 
procedures. 

Another limitation that faces the procedure is the assumption that the prepa- 
ration acts only on the system and not on the environment. This assumption 
is crucial, as we are mapping from the set of preparations on the system to the 
corresponding output states. If this assumption fails, then we would need to 
make a set of preparations that span the space of operations on the combined 
system-environment space. However, the environment can be arbitrarily large 
and we do not have any control over it. Therefore the tools given in this article 
may not be valid when the preparation affects the environment directly. When 
the preparation procedure acts on £ as well as S, the positivity of VW-map may 
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be affected. Note that, as long the effect of all preparations on £ is a constant 
for then our prescription remains valid. 

Lastly, since VW-map contains all dynamical information, we are able to 
construct B CP of Eq. Q from it. Similarly, from the correlation-memory matrix, 
we can construct B a of Eq. Q. Knowing the two we can determine B of 
Eq. (J2j|, which can be a not-completely positive map. This gives operational 
meaning to not-completely positive dynamical map as the descriptor for the 
dynamics of the system when identity preparation is made. On the other hand, 
the non-completely positive map is not experimentally determinable without 
determining VW-map. Finally, it remains an open question, when JC > 0, is B 
not completely positive? 



Methods 

The calculations in this sections are done in terms of matrix indices as M- 
map and the correlation-memory matrix K. are nontrivial tensors. We use the 
Einstein summation notation, i.e., repeated indices are summed over. Bipartite 
state of S£ is expressed with four indices with the Latin indices belong to 
S and greek indices to £. For instance, the state in Eq. ([I]) has the form 
pfa,sp = Prs ® Pa/3 + Xra,s0- A map acting on a density matrix is written 
as A rr ';ss' (pf's') = ^rr' Pr' s' A k * ss > = Prs, where A k are the Sudarshan-Kraus 
operators (see [15]). A* is the complex conjugation of A and A rs — > A sr is the 
transpose; together they give Hermitian conjugation. 



AA — map 

Let us rewrite the generalised process equation, Eq. Q, in terms of matrix 
indices 

Qrs = U re - r i a A r ' r '/- s ' s " {p r " a - s "p) U se . s ip, (14) 

where the sum over e is the trace with respect to the environment. We are inter- 
ested in the reduced dynamics of S as a function of the preparation procedures. 
Thus, we can pull the preparation map out of everything else and regard it all 
as a map acting on the preparation map: 

Qrs — \Urt;r' uP r " a ;s" s <i-_s' p\ CAr'r";«'s" ) 

J^irr' r" :ss' s" (Ar'r" -s' s" ) • (-^) 



In the last equation, the matrix M. is defined as: 

■>/3 u se,s'P 



M-rr'r";ss's" — U r£ y a P r " ct.s" Bases' 0- (16) 



Determining AA— map 

Let {p( m ) = |7r( m )) (7r( m )|} be a set of pure states that linearly span the space 
of S. There are d$ such matrices. That is, any state of S can be written as a 
linear sum of these pure states: p s = Yl m r m P ■ 
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A preparation map acting on S is a d 2 s x <P S Hcrmitian matrix. Therefore, any 

matrix in this space can be spanned by a tensor product of the basis matrices 
{pM p(n)^ which 

is a basis in for dg x d 2 s space of maps. There are dg 
elements in the basis {p( m ) <g> P(")}. We can write action of one of these basis 
element on a density operator on S as 

A {mn) {p S£ ) = \n (n) ) (tt^Ip^Itt^) (tt^| . (17) 

It is crucial to note here that {-k^ \n^) ^ 5 mn , as these vectors are eigenvectors 
of the basis elements {p(" 1 )} that do not commute. 

These preparations are can be thought of as a projection followed by a 
rotation. Action of any map on space of S acting on the SS state can be 
expresses as a linear sum 

A (p S£ ) = j2 a(mn)A(mn) ( 18 ) 

mn 

= J2 a(mn) l 7 ^) (7r (m) |p 5£ |7r (m) ) (7r (n) | (19) 

mn 

= J2°t imn) P (n) L 5£ p(™)j (20) 

mn 

_ ^ a H p N P (») 8 /IW ] (2i) 

mn 

where a*-" 1 ™- 1 are the coefficients that determine A in terms of {A^ mn ^}. p £ K m ) 
is the conditional state of the £ and p( m ) = tr[P^ m ^p 5 ] is the probability for the 
outcome p( m K 

Knowing the output states corresponding to each of these inputs, 

Q(mn) = tr£ pp(n) ^ p £|(m) {7 t] ; ( 2 2) 

along with the success probabilities p( m \ for all m,n, is enough to predict the 
output state for any preparation: 

Q = tr £ [UA(p S£ )tf] (23) 
= J2 aimn) P (m)tY £\- Up(n) ®P £](m) tf] (24) 

mn 

= ^a (mn) p {m) Q (mn) . (25) 

mn 

.M-map can be determined choosing A^ mn \ followed determining the corre- 
sponding Q( mn } and p^ m \ and standard inversion techniques [37]. Note that 
any other set of linearly independent preparation can be linearly mapped to the 



preparations given in Eq. (17), and therefore will suffice. 

Determining all Q( mn > is done by quantum state tomography. This is equiva- 
lent to carrying out dg standard QPT procedures, one each p f IM. Additionally 
measuring p( m ) is equivalent to doing quantum state tomography of p s . 
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Before moving on a simple example may be useful. For one qubit, we may 
take the following projectors as a linearly independent basis: 

p(D = I(I + (Tl ), P(2) = I(I + (T2 ), p(3) = 1(1 + a 3 ), P( 4 ) = i(l-a 1 ). 

Note that, this is a linear but not a convex decomposition: P' 5 ) = |(I + 02) = 
p(i) + p(4) _ p(2)_ The eigenvectors of pW, p( 2 \ p( 3 ), and P< 4 ) are |ar+) = 
^(|0) + |1)), \y+) = ^(|0> + \z+) = |0), and = ^(|0) - |1» 

respectively. Using these eigenvectors we can write basis elements for the maps 
that operate on the space of one qubit. For instance, 

A^\p S£ ) = \x+) (x+\ o S£ \x+) (x+\ = p(*+) \x+) (x+\ ® p £ ^ x+ \ 
A [3A) (p S£ ) = \x-) (z+\ P S£ \z+) (x-\ = V {z+) V-) (x-\ <8 p £][z+ \ 
A^ 2) (p S£ ) = \y+) {x-\p S£ \x-) (y+\ = p(*-> \y+) (y+\ ® p £ ^~\ 
and so on. 



Detecting initial correlations 

The initial state of the system is labeled by indices r" and s" . Tracing over 
everything else we can hnd the initial state of S (before preparation) from M.- 
map: 

1 s 

"■s 

This is, of course, attainable by doing state tomography at the beginning of the 
experiment, by measuring the values of p( m ) from last section. 
Next, let us the trace over the system indices r" and s" 

S r " s " M. rr i r fi - ss i s ii — Umy a p a pU ses , p = B^. ss r. (27) 

The last equation is exactly the dynamical map in the absence of initial corre- 
lations, given in Eq. Q. In other words, in the absence of initial correlations, 
QPT would yield this map. 

This means, even though the .M-map contains the information about uncor- 
rected S£ state and the correlations separately. Consider the following matrix 



composed of the matrices in Eqs. (26) and (27 1 

>-rr'r"-,ss' s" u rr' ;ss' Hr" s" 

= U re y a p' s „ s „p a/3 U* es , /3 . (28) 

The last equation is similar to the expression for the .M-map, except the state 
of the system and the state of the environment are uncorrelated. 
Writing the state of S£ in A^-map in terms of Eq. ([I]), we get 

M. r , r „. s , s „ — U rei r' a (p S „ s ,,p £ t f 3 + Xr"a;s" p)Ut e , s > p- (29) 

Now we can define the correlation-memory matrix as 

JGf-r' r" ;ss' s" J^t Tr ' r " - ss ' s " £ rr > r " - ss ' s " (^0) 
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Properties of M 
Linearity 

Mathematically, M.— map acts on the preparation map just as the dynamical 
map acts on a density operator. In fact, we are not varying the initial state 
of the system, rather the preparation procedure on that state. Therefore the 
linearity of quantum mechanics is preserved for the .M-map acting on different 
preparation procedures, i.e. 

M \ ai A {1) + a 2 A {2 A = ai MA {1) + a 2 MA {2) . (32) 

This is very much like the dynamical maps action on mixtures of states. Further- 
more, if we show that the M.— map preserves trace, Hcrmiticity, and positivity 
on its domain then all of these properties will be preserved on the state space. 
In other words for any preparation, A^ that preserves trace, Hermiticity, and 
positivity, the action of the M.— map on it will yield an output state, Q^ m \ that 
is unit-trace, Hermitian and positive. 

Trace preservation 

Let us start with the trace of M with respect to the final indices r with s: 

tT(r S )[M.\ — S rs Ai r r'r";ss's" 

= Urcj-'apf^gHfjU^g/p. (33) 

Since WU = ^ re U* es , pU re , r > a = I, then 

tr (rs) [M] = <W(Wf,? Q , s "/3 = 1 ® P S - ( 34 ) 
A preparation acting on the above matrix will yield 

tT (rs) [M}(A)=tT[A(p s )} = l. (35) 

The implication being M. preserves the trace of A{p s ). As long as the prepa- 
ration is trace a preserving operation we get a unit-trace matrix for the output 
state. 

Hermiticity preservation 

As with the case of general quantum operations, matrix M is Hermitian. This 
is easy to see by taking the complex conjugate of matrix M, 

(M rr 'r":,ss's")* = ^ U se,s' Ps»/3 ,r" a U re,r' a, 
a/3e 

= ■M ss 's";rr"r' ■ (36) 

The complex conjugate of M is not only the transpose of M, but each element 
of M is also transposed. Hence M is a Hermitian matrix. 
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Positivity of A4-map 

The .M-map is composed of a unitary matrix operating on a density matrix. 
Then we can take the square root of the density matrix to get 

KA , „ , „ TJ , a/ nS£ ,/ n S£ rr* 

jvi r r'r";ss's" — ^re,r'aV P r"a,a~/VI J a^.s" fi u se,s' f! 

= M^, r „M^„. s , (37) 

where M = U \J p S£ and /i = aje. We have written the A^-map in operator sum 
representation, hence it is completely positive. Where are the Sudarshan- 
Kraus operators |15[ I16j. This means, the A^-map acting on any preparation 
procedure will lead to a physical state. This was not the case when a standard 
QPT procedure is carried out on initially correlated S£ states. The action of 
.M-map can now be written as 

M(i) = ^MMM" t . (38) 

The properties shown above are precisely the conditions for a generic quan- 
tum operation to preserve trace, Hermiticity and positivity. Therefore AA— map 
preserves the attributes on the preparations, which in return will preserve these 
attributes on the states. 
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